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We present a new approach for obtaining very precise integration results for infrared vertex and 
box diagrams, where the integration is carried out directly without performing any analytic in- 
tegration of Feynman parameters. Using an appropriate numerical integration routine with an 
extrapolation method, together with a multi-precision library, we have obtained integration re- 
sults which agree with the analytic results to 10 digits even for such a very small photon mass as 
10~ 150 GeV in the infrared vertex diagram. 
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1. Introduction 

In the field of the particle physics more and more precise computation is required for the 
progress of the accuracy of the future collider experiments. For this reason, the computation of 
one- and higher loop corrections is mandatory. 

General solution has been known for one-loop integrals [[Ij] and there are several excellent pro- 
grams developed in this decade such as FF[|2|], FormCalc-LoopToolsQ, XLOOPS-GinacQ. 
In parallel, there have been continual efforts to carry out loop integrals by fully or semi-fully nu- 
merical methods. In early 1990 Fujimoto et al. [jD, |6[ f7|, ||] have proposed several excellent methods 
and have shown numerical results. In early 2000 Kurihara et al.^\ have developed the numerical 
contour method and have shown numerical results. 

Since 2003 we have been developing a method which enables us to carry out the loop integrals 
in a completely numerical way. In this numerical method, to prevent the integral diverging, we put 
ie in the denominator of the integrand. Here £ is a real positive constant and it is not necessary to 
be infinitesimal. Thanks to this ie, the denominator does not vanish and we get the numerical result 
of 1(e) for a given e. Calculating a sequence of /(£/) varying £/ = £o * (const. )~ l (I = 0,1,2, • • •) 
and extrapolating them, we can get the final result of the integration in the limit of e — > 0. So far in 
this method we have calculated several one-loop and two-loop integrals and reported the numerical 
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Figure 1: One-loop vertex diagram 



results in [[H], 11, 12, 13]. In these demonstrations it has been clearly shown that the extrapolation 
method is very efficient for not only one-loop but also two-loop diagrams. 

In the numerical calculation of the loop integrals, it is also very important to confirm whether 



the method can handle the infrared singularities. In [12], we have calculated the loop integrals 
with the infrared singularities with a photon mass of up to 10~ 15 GeV. In the calculation we used 
the double precision arithmetic for the photon mass of up to 10 -6 GeV. However for the photon 
mass of below 10 GeV the double precision is not enough and we use the quadruple precision 
arithmetic. In the quadruple precision calculation the relative error is given the order of 10~ 7 . 

In this paper, we propose a new approach for our numerical method. In this approach we 
include a precision control technique in addition to the extrapolation method for the loop integral. 
We show that a precision control is mandatory to get the stable results of one-loop diagrams which 
have infrared singularities with smaller photon mass below 10~ 30 GeV. 

The layout of this paper is as follows. In § ^ we give formulae of the loop integrals we consider 
in this paper. In § ||| we give an explanation of a new approach of our numerical method. Numerical 
results are shown in § ||. Finally in § || we will summary this paper. 

2. Loop Integrals 

2.1 One-Loop Vertex 

For one-loop vertex diagram in Fig. [TJ the loop integral we consider in this paper is 



/ dx dy-, (2.1) 
Jo Jo D 



where 

\2 2 



D= -xys+(x + y) m 1 + {\ - x-y)X z . (2.2) 

Here s denotes squared central mass energy and m and A are a mass of external particles and 
a fictitious photon mass respectively. We introduce A so as to regularize the infrared singularities. 
Replacing s by s + is in ( 2.2), the real part of the integral ( |Q| ) becomes 
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f D 

9v = / w5 . - ~ dxdy, (2.3) 

JsD 2 + e 2 x 2 y 2 

and the imaginary part becomes 

z= Lw^h? dxdy - (2A) 

Here 5 is the triangular region as < x, < y nn&x+y < 1. Corresponding analytic formulae 
for one-loop vertex diagrams are given in ^ ||, g]. 

2.2 One-Loop Box 

For one-loop box diagram in Fig. ||, the loop integral we consider in this paper is 

1 rl-x rl-x-y l 

dx dy dz-=~, (2.5) 

o Jo Jo O z 

where 

D = -xys-tz{\ -x-y-z) + {x + y)X 2 + {\ -x-y-z){\ -x-y)m 2 +z(l -x-y)M 2 . (2.6) 
In ( ^ ) s denotes squared central mass energy and m and M are external masses of particles. 




Figure 2: One-loop box diagram 



When we replace s by s + ie in ( |23| ), the real part of the integral becomes 

91= / m , 2 2 2 , 2 dxdydz, (2.7) 
and the imaginary part becomes 

f lexyD 

3 = / — ^ 7 ». n dxdydz. (2.8) 

Js (D 2 + e 2 x 2 y 2 ) 2 

Here 5 is the tetrahedron region as < x, < < z and x + y + z < 1 - Corresponding analytic 
formulae for one-loop box diagrams are given in [|6|]. When we put M = }X, the scalar integral for 
one-loop box diagram becomes 
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analy " C s ^/[-t + n 2 - m 2 ) 2 - 4m 2 t 



x In 



3. Numerical Method 



{-t + H 2 + m 2 + ^/(-t + n 2 - m 2 ) 2 - 4m 2 t) 2 
V X 2 ) \m 2 ]X 2 ■ { ' 



The basic idea of our numerical method is the combination of an efficient multi-dimensional 
integration routine and an extrapolation method. Using this method we have shown several nu- 
merical results of loop integrals and that they show good agreements with analytical results. In 
the calculation so far we have used double precision arithmetic or quadruple precision arithmetic 
if necessary for one-loop and two-loop integrals. However for the infrared divergent diagram even 
with the quadruple precision arithmetic it becomes harder to get a result with an enough accuracy as 



the value of the A becomes much smaller such as 10~ 15 GeV[12]. Therefore in our new approach 
we include the concept of a precision control. 

In the following subsections, we give a very simple explanation of the multi-dimensional inte- 
gration routine, an extrapolation method and a multi-precision library we used in our new approach 
respectively. 

3.1 Multi-dimensional Integration Routine 

We use DQAGE routine in the multi-dimensional integration. It is included in the package 
QUADPACK and it is a globally adaptive integration routine. Though it is a routine for one- 



dimensional integration, we use it iteratedly Q15Q . 



3.2 Extrapolation Method 

Although there are several ways to accelerate the conversion of the sequences. Of them, the 



e algorithm by P.Wynn[16] is recommended as the best all-purpose method for slowly converging 



sequences. It is an implementation of Shanks' transformation [ 17] and is represented in a recursive 
formulae. 

3.3 Multi-Precision Library 

When we know correctly whether the precision arithmetic used in the calculation is sufficient 
or not, we can select the necessary precision arithmetic and the results become reliable. J.Fujimoto 



et. al have discussed the importance of a precision control in Q18p. We use HMlib[]19[] as the 
multi-precision library in our new approach because it gives an information of the lost-bits during 
the calculation and we can guarantee the precision of the results. Taking the calculation of one- 
loop vertex diagram as an example, we show the information of the lost-bits supplied by HMlib in 
Table [!]. 
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Table 1: In P-precision presentation implemented in HM1 ib, a sign bit is 1 bit and an exponent bit is 15 bits 
and a mantissa is (32 x P — 16) bits. When P — 4 (quadruple precision), the mantissa becomes 1 12 bits. In 
the Table when A is 10~ 23 GeV, the maximum number of lost-bits is the same as the number of bits of the 
mantissa. 



X [GeV] 


Average lost-bits 


Maximum lost-bits 


10 -20 


88 


92 


io- 21 


98 


102 


10 -22 


108 


112 



Table 2: Numerical results of one-loop vertex diagram with y/s = 500 GeV, t = — 150 2 GeV 2 , m = m e = 
0.5 x 10~ 3 GeV. The upper is the result of Real part and the lower is one of Imaginary part. 



A [GeV] 


Numerical Results 


Precision 


Analytical Results [quadruple precision] 


Agreement 


io- 3U 


-0.150899286980769753D-01 ± 0.771D-26 
0.189229839615898822D-02 ± 0.124D-25 


8 
8 


-0. 1 5089928698048229 151911 707856807798E-0 1 
0. 1 892298396155253893453366778083205 1E-02 


12 
12 


10 -60 


-0.30359395256285495 1D-01 ± 0.178D-16 
0.3628406655 14227622D-02 ± 0.896D-13 


16 
16 


-0.30359395256226015422329307055958645E-01 
0.36284066551349654484509779501460548E-02 


12 
11 


10 -8U 


-0.405390396284235075D-01 ± 0.580D-15 
0.478581216125478532D-02 ± 0.401D-12 


16 
16 


-0.40539039628344539602607706522058889E-01 
0.4785812161 1214398184493853981879066E-02 


11 
10 


10 -uo 


-0.608983283726997427D-01 ± 0.556D-15 
0.710062317325786663D-02 ± 0.415D-12 


32 
32 


-0.60898328372581587963164505454259761E-01 
0.7 100623 1 7309438855 844620029427 1 6 1 00E-02 


11 
10 


10 -uo 


-0.761677949309069452D-01 ± 0.931D-15 
0.883673 1431 858014 14D-02 ± 0.260D-13 


32 
32 


-0.76 1 67794930759374233582 1 045 1 6222 1 89E-0 1 
0. 883673 1 432074 1 00 1 1 34438 1 1 4507364254E-02 


11 
10 


10 -1M 


-0.812576170752810666D-01 ± 0.549D-10 
0.941543418501223556D-02 ± 0.109D-11 


32 
32 


-0.81257618347269926993757808016127325E-01 
0.94154343249672244271406905030553399E-02 


7 
7 



4. Numerical Results 

4.1 One-Loop Vertex 

The numerical results are shown in Table |2[ Results are compared to analytic results evaluated 
by the formulae in [Jsj] . As an example, in Table || we show the parameters we used in the cal- 
culation which are key and limit used in DQAGE, starting £ and ending £ for an extrapolation 
method. 

4.2 One-Loop Box 

The numerical results of one-loop box diagram shown are shown in Table |j. The results are 



compared to the analytic results evaluated by ( 2.9) in a quadruple precision. As an example, in 



Table |5| we show the parameters we used in the calculation which are key and limit used in 
DQAGE, starting e and ending £ for an extrapolation method. 

4.3 Elapsed Time 

In Table |6] and ^, the elapsed times required in the several typical calculations are shown. 
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Table 3: Parameters used in DQAGE and in an extrapolation method for the one-loop vertex diagram with 
sfs = 500 GeV, m = m e = 0.5x 10~ 3 GeV, A = 10~ 150 GeV. The upper is the parameters for the Real part 
and the lower is ones for Imaginary part. 



key x 


limit x 


£.Sf£!rt 


£end 


key y 


limity 






1 


600 


0.36572620D+04 


0.28485158D+03 


1 


600 






1 


600 


0.36572620D+04 


0.28485158D+03 


1 


600 







Table 4: Numerical results of the loop integral of one-loop box diagram with ^/s = 500 GeV, t = -150 2 GeV 2 , 
m = 0.5 x 10~ 5 GeV, M = 150 GeV. This is results of Real part. 



A [GeV] 


Numerical Results 


Precision 


Numerical Results [quadruple precision] 


Agreement 


i<r 14 


-0. 1927861 102670278D-06 ± 0.314D-14 


4,4,2 


-0. 1 92786 1 1 2243964 1 1 60689577 1 777 1 95708E-06 


8 


10 -a> 


-0.2472486348234972D-06 ± 0.586D-15 


4,4,2 


-0.24724863525991758999379683385554222E-06 


9 


io- 25 


-0.30171 1 1253761463D-06 ± 0.11 1D-13 


4,4,2 


-0.30 1 7 1 1 1 5 8275 87 1 0639 1 8635949939 1 3029E-06 


7 


i(r*> 


-0.3562028882831722D-06 ± 0.867D-10 


4,4,2 


-0.356 17368 129 1 82453784347506602271 837E-06 


4 



Table 5: Parameters used in DQAGE and in an extrapolation method for the one-loop box diagram with 
sfs = 500 GeV, t = -150 2 GeV 2 , m = 0.5x 10~ 3 GeV, M = 150 GeV and A = 10~ 25 GeV. 



key x 


limit x 






key y 


limity 






key z 


limit z 






1 


400 


0.36572620D+04 


0.28485158D+03 


1 


400 






1 


100 







Table 6: Elapsed time required in the calculation of the one-loop vertex diagram 



A [GeV] 


Real Part 


Imaginary Part 


Precision 


CPU 


10 -30 


1.8days 


1.8days 


8,8 


Xeon 3.06GHz 


1Q -150 


6.8days 


6.7days 


32,32 


Xeon 3.06GHz 



Table 7: Elapsed time required in the calculation of the one-loop box diagram 



A [GeV] 


Real Part 


Precision 


CPU 


io- 15 


1.5days 


4,4,2 


Opteron 2.2GHz 


lO" 25 


3.0days 


4,4,2 


Opteron 2.2GHz 
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5. Summary 

We have shown that a new numerical approach presented in this paper is applicable to the 
scalar one-loop vertex and box diagram with infrared singularities. The new approach consists of 
an integration routine, an extrapolation method and a precision control. Several numerical results 
are shown and they agree with analytic ones. All the calculation of the one-loop vertex diagram 
are carried out in at least quadruple precision arithmetic and some are in 8-, 16- and 32-precision 
arithmetic. From these demonstrations, it is confirmed that a precision control plays an important 
role in the calculation of the loop integrals with the infrared singularities. Although higher precision 
arithmetic costs a lot of CPU time, it supplies the information of the lost-bits occurred during the 
calculation and this is extremely valuable in getting the reliable results. To reduce CPU time we 
will be able to use the parallel computing technique. 
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